arXiv: 1505.05517v2 [hep-th] 5 Aug 2015 


P — v Criticality in Quasitopological Gravity 


Robie A. Hennigarj^] W. G. Brennaj)] and Robert B. Manrj^] 

Department of Physics and Astronomy, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada 

We investigate the thermodynamic behaviour of AdS quasitopological black hole solutions in the 
context of extended thermodynamic phase space, in which the cosmological constant induces a pres¬ 
sure with a conjugate volume. We find that the third order exact quasitopological solution exhibits 
features consistent with the third order Lovelock solutions for positive quasitopological coupling, 
including multiple reentrant phase transitions and isolated critical points. For negative coupling we 
find the first instances of both reentrant phase transitions and thermodynamic singularities in five 
dimensions, along with other modified thermodynamic behaviour compared to Einstein-AdS-Gauss 
Bonnet gravity. 


I. INTRODUCTION 

Black hole thermodynamics has maintained its pop¬ 
ularity, partly because black holes make for compelling 
and elegant thermodynamical systems whose behaviour 
can be linked to the nature of quantum gravity. Of partic¬ 
ular relevance to this work are the asymptotically Anti 
de-Sitter (AdS) black holes, which have significance in 
various proposed gauge-gravity dualities. 

Extended phase space thermodynamics, in which the 
cosmological constant A is regarded as a thermodynamic 
variable [T| analogous to pressure EHH], has recently be¬ 
come of considerable interest as the proper venue for a 
complete thermodynamic description of AdS black holes 
m- The mass of an AdS black hole is then understood 
as the enthalpy of the spacetime [3]. It was soon re¬ 
alized that the critical exponents of the 4-dimensional 
Reissner-Nordstrom AdS black hole are the same as those 
in the van der Waals liquid-gas phase transition 12] . 
completing an analogy explored previously in a more re¬ 
stricted context [2QH23]- Subsequently a broad range of 
new thermodynamic behaviour was discovered, including 
reentrant phase transitions, tricritical points, and more 
general van der Waals behaviour with standard critical 
exponents mmm , with an extensive review of these 
issues in the context of rotating black holes presented in 

m- 

One of the most recent discoveries of interest was ob¬ 
tained in 3rd-order Lovelock gravity [25; :H]. Multiple- 
reentrant-phase-transition behaviour as well an isolated 
critical point with non-standard critical exponents for a 
D = 7 third-order (cubic) Lovelock asymptotically AdS 
black hole were found. This latter phenomenon is the 
first example of a critical point with non-standard criti¬ 
cal exponents obtained in a geometric theory of gravity. 
It was subsequently shown that it occurs for topological 
black holes in any I\ = (2k + l)th-order Lovelock gravity 
theory, provided the Lovelock couplings are appropriately 
adjusted [Hj. This corresponds to a double swallowtail 
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behaviour in the Gibbs free energy, giving rise to two 
first-order transitions between the small and large black 
hole phases. 

Our key interest for this paper is to see if these more 
recently discovered phenomena can occur in a broader 
context. To this end we explore the extended phase space 
thermodynamics of black holes in cubic quasitopological 
gravity. Quasitopological gravity was proposed pfoT - FTTl as 
an extension to the Lovelock [35] higher curvature terms. 
While the Lovelock terms are topological invariants be¬ 
low a certain critical dimensionality, the quasitopologi¬ 
cal terms are still active in some lower dimensions. The 
trade-off is that when spherical symmetry does not hold, 
the quasitopological terms produce higher than second 
order derivative terms in the Einstein field equations. 

With cubic quasitopological gravity we can obtain field 
equations in D = 5 similar to those of cubic Lovelock 
gravity. We will examine whether we can likewise obtain 
similar thermodynamics but in fewer dimensions. One 
particularly compelling notion for a D = 5 black hole 
is that the basic idea behind AdS/CFT will link it to 
a D = 4 field theory, and so perhaps universality rela¬ 
tionships with well understood D = 4 thermodynamic 
systems will become apparent. We find that indeed new 
features are present in D = 5, particularly the first in¬ 
stance of reentrant phase transitions, a phenomenon so 
far seen only in D > 6 [27 . 

We consider quasitopological additions corresponding 
to 3 rd -order curvature corrections to Einstein-Maxwell- 
Gauss-Bonnet gravity that maintain second-order field 
equations with respect to the metric under conditions of 
spherical symmetry. We employ the action 

' = fd D *V=i (-2A + c, + (d _ 3) a (d _ 4) £2 

+ (D-3KD-6) £ ^^ F “") (1) 

where the couplings have been redefined from their orig¬ 
inal formulation [451447] to allow a direct and straight¬ 
forward comparison with Lovelock gravity. Here D is 
the number of dimensions, F jiV = d^A^, /r and A are 
the correction terms’ coefficients, C\ = R is the Ricci 
scalar, £ 2 = R^ vl sR^ lS — 4 + R 2 is the Gauss- 
Bonnet Lagrangian, and £3 is the 3 rd order quasitopo- 
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logical gravity term. This term has the form 


II. EXACT SOLUTION 


£3 = 


2D-3 


3D 2 - 15 D + 16 


R„\ p Rv T B a RS 


+ 7 -- 7 ( ^ R l , a ^R liai ' ri R 

(2D — 3)(D — 4) V 8 ^ d 

- 3 (D - 2)R liai/ pR^\RP T + 3D ■ R^R^R^ 

+6(D - 2 )R^R a u R^ 


3(3£>~4) , 3^.3 

tX/j, il a IX “h c -tl 
Z o 


( 2 ) 


and is only effective in higher dimensions (U > 4), while 
becoming a surface term in 6 dimensions 145| . 

Of particular importance to this work is the scaling 
of the higher-curvature coupling coefficients. Written in 
this form, A and /i are dimensionful, since the Riemann 
terms obey scaling ~ L~ 2 , where L is some lengthscale 
set by the metric coordinates. Therefore, A ~ L 2 and 
fj, ~ L 4 . Later, when we use an Eulerian argument to 
construct an extended Smarr formula from a 1 st law, the 
scaling of these terms will set the coefficients in the Smarr 
formula. 


The metric ansatz we employ is 

ds2 = ^f( r ) dt2 + ^2g^ +r2dfl l ( 3 ) 

where dfl 2 is the metric of a constant curvature hyper¬ 
surface 



I d0 2 2 + n sin2 9 i de i j 

\ i =3 3=2 ) 


(4) 


The parameter k = —1,0,1, corresponds to hyperbolic, 
fiat, and spherical geometries, respectively. Topological 
black holes are obtained by putting appropriate identi¬ 
fications on the former two [ 501 . We parameterize the 
Maxwell gauge field with a function h(r) as 

A = q^rh(r). (5) 

Substituting this ansatz into the action, using 


In the remainder of this work, we will first generalize 
the exact solution found in [49] to both D = 5 as well as 
D = 7. We will compare the resulting thermodynamic 
behaviour of the D = 7 solution with previous work from 
cubic Lovelock gravity |28| . and we will discuss any new 
features that arise as we take D = 5. 


(D - 1 )(D - 2 ) 

2L 2 


( 6 ) 


as the cosmological constant and simplifying using inte¬ 
gration by parts to remove terms proportional to deriva¬ 
tives of f(r), we find 


I = 


1 

Z# 


d D ~ l x 


dr \fl{{ iD - 2yD ~' 





+ 2 q 2 r 


D-2 9_ 

f 



(7) 


where a prime denotes differentiation with respect to r 
and k = g(r) — L 2 k/r 2 . Note that in eq. (|7| the deter¬ 
minant of the angular part of the metric has been sup¬ 
pressed, since it is unimportant in obtaining the field 
equations. 

We substitute f(r) = N 2 (r)g(r ) and vary the action 
with respect to g(r), N(r) and h(r) to obtain the field 
equations 


(D — 2)r D ~ 1 



^1 — k + 



(8) 

^2 d 3^1 
L 2 L 4 ) 

/ 


2 a 2 r D ~ 2 

V [(<*)' 

2 

(9) 

(r D ~ 2 (rh)'\ 
l N ) 


(10) 


Clearly, N = 1 solves equation ([8]) and this solution 
will be taken for the remainder of the paper. Finding 


h(r) from equation (10) and solving for the gauge field 
from © , we find 


At — 


I (D-2) q 
2 (D - 3) r D ~ 3 


( 11 ) 


Solving the remaining differential equation, we obtain the 
fundamental relationship 


1-8 + l ? k2 - = T 


L 2 m q 2 L 2 

D-l — r 2{D-2) 


( 12 ) 


where 


L 2 

R = -v (F(r) - k) 


(13) 


and F(r) = r 2 g(r)/L 2 is defined such that the metric 
reads: — F(r)dt 2 + ■ ■ •. Note that if we were to use the 
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function g(r ) as defined previously, we would be exclud¬ 
ing factors of L from the mass and charge. Since in ex¬ 
tended phase space we consider L to be a parameter that 
is varied, it is more convenient to write expressions in 
terms of F(r). Note also that the mass parameter m and 
the charge parameter q have been defined such that we 
recover the RN-AdS solution in the form presented in [T3j 
when we take /i = A = 0 . 


III. SMARR RELATION AND 1 st LAW 


The remaining potentials can be calculated by requir¬ 
ing consistency of the first law. Defining the pressure 
as 


A _ (D- 1)(D — 2) 
87t 167tL 2 


a straightforward calculation yields 


V = 


Wfl-2 D — l 

D- I rfl 


for the thermodynamic volume and 


( 20 ) 


( 21 ) 


When considering the thermodynamics of these black 
holes it is useful to work explicitly with eq. (12). The 
utility of using this expression is that one is not restricted 
to considering only one branch of the solution. As a result 
of this, it is natural to expect that all thermodynamic 
quantities should limit to the standard D —dimensional 
RN-AdS results as g = A = 0. 

We now move on to considering the thermodynamics 
for these black holes in extended phase space, where all 
dimensionful parameters are considered to be thermody¬ 
namic variables EH- We employ the extended first law, 


= up- 2 {D-2) d _ 5 f _ 8t rr h T \ 

16tt h \ D - 4 ) ’ 

_ uj d _ 2 (D- 2) 2 n_ 7 ( 12nr h T \ 
16tt h { D - 6 ) 


( 22 ) 

(23) 


for the potentials conjugate to the couplings. One finds 
that these thermodynamic quantities satisfy the Smarr 
relation 


(D -3 )M=(D- 2 )TS - 2 PV + (D - 3) 
+ 2A4- a + 4/iT /; 


dM = TdS + VdP + + T A dA + (14) 


where 


,, (D- 2) 

M = ——— -UJ D _2 m 

l07T 

„ V / 2( j D-2)(D-3) 

Q = -o- UD-2 q 


Sn 


as in Enn], where 


ujd = 




(15) 

(16) 

(17) 


is the surface area of a D-dimensional sphere. 

The potentials introduced by the quasitopological and 
Gauss-Bonnet parameters follow from the procedure in¬ 
troduced in [2S]. The entropy of this solution was found 
to be 


r^~ 2 uj D -2 ( 2 (D-2)Xk 3(D-2)gk 2 \ 

4 ^ + p-4)r 2 + (D-6)ri ) 

(iS) 

and the temperature can be found from 


T = 


/ dF{r)/dr\ 
V 4t t ) 


r= 


Th 


which is consistent with Eulerian scaling. 

We are now in a position to examine these black hole 
solutions for critical behaviour in extended phase space. 
In the following sections we will adopt an approach which 
is similar to that used in 128j insofar as we will employ 
dimensionless thermodynamic parameters to simplify our 
investigation. However (unlike the analysis performed in 
[28] ) we will not restrict ourselves to positive values of 
the coupling [i but will instead explore the entire param¬ 
eter space that is consistent with the entropy, pressure, 
and stability constraints (or in other words a ghost-free 
vacuum) discussed in the following sections. The case of 
positive quasitopological coupling will be considered in 
section |IV[ while section [V] will focus on negative qua¬ 
sitopological coupling. 


IV. POSITIVE \x THERMODYNAMICS 

Analogous to [281. we define the dimensionless thermo¬ 
dynamic quantities^ 


Th = vfi 1/4 , 


q = 


ql(D-2)(D-3) _ (D _ 


(D- 3)/4 


t = (D — 2 )h 1 /4 T. 


(24) 


where is the outermost root of F(r) = 0. In practice, 
to compute dF(r)/dr it is easiest to differentiate eq. (12) 
and solve for the derivative term. Doing so yields 


T = 


1 


47T (r 


2 k\r- 2 h 


3 k 2 g) 


(d - 

L 2 


+ {D- 3 )kr\ 


~\~{D — 5)k 2 Xrh + 


(D - 7)k 3 g (D - 3 )Q 2 


rh 


r 2D-9 
' h 


(19) 


where we have performed a rescaling of the charge to 
make comparisons with [28] more convenient. Note that 


1 Note that the definition of the dimensionless charge differs from 
that in m , while all other dimensionless quantities are the same. 
This adjustment was made to the charge to account for the dif¬ 
ferences in coupling to the electromagnetic field. 
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from (21) rh will be non-linearly related to V. It is more 
convenient to work with the specific volume v, than V. 
Qualitatively this has no effect on the results we obtain, 
though care must be taken upon integration over v (for 
example in computations involving the equal-area law) 
since the measure of integration will contribute. 

In terms of these quantities, along with the parameters 


a = ~7= P = I * * 4 VpP, 
\/M 


(25) 


we can write the dimensionless equation of state in the 
form 


_t (D-2)(D-3)k 2 kta (D — 2)(D — 5)k 2 a 

P v 4nv 2 v 3 4nv 4 


+ 3 kH (D-2)(D-7)k q 


47TU 6 * 


2(D-2) ' 


(26) 


We see that eq. (26) is identical to the equation of state 
found in the study of 3 rd order Lovelock gravity (c.f. 
eq. (4.7) from JB])- Hence the same thermodynamic 
behaviour seen previously for D > 7 - particularly mul¬ 
tiple re-entrant phase transitions [25] and isolated critical 
points [44] - will be seen here as well. 

In this case we are permitted to study the D = 5 be¬ 
haviour, since quasitopological gravity is well-defined in 
five spacetime dimensions. For D = 5, some terms in 
( |26[ > vanish and others reverse in sign, and so we must 
directly analyze this particular case. We shall carry out 
this analysis in the next subsection. 

Next we compute the Gibbs free energy, G = M — TS 
for this solution (in dimensionless form): 


9 = 




(3-D)/4 


W_D-2 


-G. 


(27) 


Explicitly in the case k = ±1, this is given by 


9 = 


167 T (v 4 + 2 akv 2 + 3) 


3 (D - 2)k 


,D-7 


3 (D - 2 )(D - 8 )av D ~ 5 


(D — 6 ) 


a(D — 8 )i> 


D — l 


607 TpV 


D—l 


D-4 


{D- 1){D- 6 ) 


4 (D + 3 )kv D ~ 3 


+ 


4(D — 3) (v 4 + 2 akv 2 + 3) v D 3 


- kv D+1 + 


(D — 6) (D — 4) D- 6 

24-Kakpv D+1 4irpv D+3 


2 (D - 2 )a 2 kv D ~ 3 
D- 4 


(£>-!)(£>-4) (D — 1)(U — 2)J 


(2D - 5)v 4 2(2 D - 7)akv 2 3(2 D - 9) 


D-2 


D-4 


(D- 6) J ' 


(28) 


Direct comparison between eq. (28) and eq. (4.9) in [25] 
reveals that the two expressions are identical. That is, 
the dimensionless Gibbs energy in quasitopological grav¬ 
ity is the same as the dimensionless Gibbs energy for 3 rd 
order Lovelock gravity. 

We see from eq. (18 1 that for certain choices of pa¬ 


rameters the entropy can be negative. In the thermody¬ 
namic considerations that follow, we consider only black 
holes with positive entropy to be physical. For this rea¬ 
son, it is important to clarify the restrictions that deter¬ 
mine the positivity of entropy. The result is trivial in 
the k = 0 case, as the entropy is always positive. In the 
cases k = ±1 the condition will depend on the number 
of spacetime dimensions. Substituting the dimensionless 
volume v from eq. (24) into the expression (18) for the 


entropy and solving for the zeros gives 


v± = 


(D-2 

\J 


—ka ± 


3 (£>- 4)2 
(D-2)(D-6) 

(29) 

for the roots with v > 0. For D > 6 , the result is the 
same as for 3 rd order Lovelock gravity [28]: for ka > 0 
the entropy is always positive, whereas for ka < 0 , the 

entropy is positive provided |a| < ^ J ■ -* 41 
case D = 5 positivity of entropy is satisfied for v > v+, 


I 

regardless of the sign of ka. 

The requirement that asymptotically AdS regions ex¬ 

ist for all branches of the solution was discussed in [28] . 

We examine the asymptotic behaviour of the discrimi¬ 
nant of eq. (12). If the discriminant is asymptotically 


(at large r) greater than zero, all three possible branches 
have real-valued solutions. However, if the discriminant 
is less than zero, only one of the branches will take real 

values at large r. This corresponds to the case where 
only the branch that does not reduce to either of the 
Gauss-Bonnet solutions as a 3 —> 0 has AdS asymptotics. 

Taking the asymptotic limit, the equation defining k 
becomes 


M^oo) — 1 ^00 + £ 2^00 L 4 ^ 00 — ^ (30) 


and we find that its discriminant yields the inequality 
18/rA 4A 3 A 2 4p 27p 2 


L 6 


L 6 L 4 L 4 


L 8 


> 0 


(31) 


which is quadratic in 1 /L 2 after multiplying through by 
L 4 . Writing L 2 = (D — 1 )(D — 2)/ 167 rP and substitut¬ 
ing the dimensionless parameters p and a in eq. (25) we 
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v_ - V+ 



FIG. 1. A plot of the regions bounded by p+ and p~ in 5 
dimensions. Outside of the region, in the majority of the plot, 
we have only one branch which has a valid AdS limit, while 
inside the region all three branches have valid AdS limits. 


obtain a discriminant of 


A = 


7r 2 p 2 


(D — 2) 2 (D — l) 2 

7r 3 p 3 

(D — 2) 3 (D — l) 3 
6912?rV 


(16a 2 - 64) 
(1152a- 256a 3 ) 


(d — 2) 4 (n — i) 4 

which equals zero when 
(D-l)(D-2) 


P± = 


1087T 


9a - 2a 3 ± 2 (a 2 - 3) 


3/2' 


(32) 


(33) 


denoting a transition in the number of solutions with AdS 
limits. This is exactly the same as for the cubic Lovelock 
case, but with the important distinction that we can now 
examine the behaviour in 5 dimensions, as seen in Figure 

□ 

We can examine the stability of the AdS vacuum to 
show that all regions of Figure [T] possess at least one 
branch with a ghost-free AdS vacuum. We do this by 
applying the condition 515] 


h{h oo) — 1 R a o + j^^oo (34) 

on the ghosty graviton, which will yield h'(Roo) > 0 if 
the graviton is a ghost. This is because h'(t too) appears 
as the prefactor to the standard equations of motion for 
perturbation about the AdS vacuum, and determines the 


sign of the kinetic term in the propagator for the gravi¬ 
ton. Substituting R as a solution of 


0 — 1 — Rqo + 


^ ~2 
L 2 °° 


JLR 3 

L 4 00 


which is the condition for AdS asymptotics, we find that 
h' (Roc) < 0 for two of the three branches when both the 
conditions p > p^ and p < p + are satisfied (i.e. this is 
the region where three AdS limits exist). Furthermore, in 
the region with only one valid AdS limit, h'(Roc) < 0 and 
so we again have a solution with a ghost-free vacuum. 

Another important consideration for the thermody¬ 
namics of these black holes is the existence of a ther¬ 
modynamic singularity. It was shown in [28] that such a 
singularity exists when k = —1 and when 


dp 

dt 



= 0 


(35) 


for D = 7 Lovelock gravity. In view of the preceding anal¬ 
ysis, the same situation holds in quasitopological gravity. 
Evaluating the condition above for the case k = ±1, one 
finds that the real solution for the specific volume given 

by 


v s = \J—kot ± \Ja 2 — 3 (36) 

results in a thermodynamic singularity. It is clear that 
only for ka < 0 can a thermodynamic singularity occur. 
At this value of the specific volume the temperature di¬ 
verges. However, one can find a pressure for which the 
temperature is in fact well behaved at this singularity. 
In [21] this was done by solving 


dp 

dv 




= 0 


(37) 


for the temperature, inserting it and v s into the equation 
of state to obtain the pressure at this singular point, p s . 

Since the result here is identical to that found in the 
case of 3 rd order Lovelock gravity, we see that cubic qua¬ 
sitopological gravity has the same isolated critical points 
as 3 rd -order Lovelock gravity. These occur for a = y/3 
(ensuring their entropy is positive), and correspond to the 
existence of two swallowtails, giving rise to two first-order 
phase transitions between small and large black holes. 
This phase transition is characterized by non-standard 
critical exponents in the phase diagram [44]. 


A. Thermodynamics in D = 5 

In the previous section, we have shown that the ex¬ 
tended phase space thermodynamics in 3 rd order qua¬ 
sitopological gravity is identical to the thermodynamics 
of 3 rd order Lovelock gravity, as considered in [28]. For 
this reason, we now study the results in D = 5 to see if 
there exist any new thermodynamic features not present 
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in the Lovelock case, and to see which of the thermody¬ 
namic results from Lovelock gravity are consistent with 
the constraints of D = 5. 

The first thing we note is that thermodynamic singu¬ 
larities in D = 5 are not physical. They correspond to 
black holes with negative entropy, since upon comparison 
of eq. (29) with eq. (36) we find v + > v s . Hence cubic 
quasitopological gravity has no special isolated critical 
points in five dimensions. 

Moving on to a more detailed study, we note that due 
to the form of the equation of state, there is no interesting 
thermodynamic behaviour when k = 0. Henceforth we 
examine only spherical and hyperbolic black holes. 


1. Spherical 


In the case k = 1, the conditions for critical points, 


d P = 0. *1=0. 

dv ’ dv 2 


(38) 


can be solved analytically, although the resulting expres¬ 
sions for t c ,v c and p c are not illuminating. The result is 
that there exists one critical point which satisfies the con¬ 
straint conditions provided a > ~ 1.9909 (smaller values 
of a result in a p c that is smaller than p_). 

We can calculate the critical exponents corresponding 
to this critical point by expanding the equation of state 
in terms of the parameters 


w = 


- L 



(39) 


In terms of these parameters our expansion takes the 
form, 

— = 1 + At + Btoj + Coj 3 + 0(toj 2 , w 4 ) (40) 

Pc 


where A, B and C are complicated a-dependent coeffi¬ 
cients. However, the important feature of this expansion 
is the form, which is identical to the cubic Lovelock case, 
and to the van der Waals fluid. Because the equation of 
state entirely governs the macroscopic behaviour of the 
system, the critical exponents are therefore given by the 
standard mean field theory results: 

a = 0, 7 = 1, = 3. (41) 


2. Hyperbolic 

We begin by mentioning that there are no isolated criti¬ 
cal points found in D = 5. In the case of Lovelock gravity, 
isolated critical points were found in D = 7 and higher 
and correspond to the case where the critical point occurs 
at the thermodynamic singularity. However, as we dis¬ 
cussed in the previous section, when D = 5 and k = — 1 


the thermodynamic singularity occurs when the entropy 
is negative. Thus the black holes exhibiting the isolated 
critical point are unphysical. 

In fact, the D = 5 hyperbolic black holes appear to 
exhibit no critical behaviour at all. When one solves for 
v c it is found that there are solutions for two critical 
points. However, when one compares the solution for v c 
and p c against the constraints on pressure and entropy, 
it is found that the solutions violate either one or both 
of these constraints. Thus, the critical behaviour for the 
hyperbolic black holes is unphysical. 

One can consider the thermodynamics outside the 
pressure wedge of Figure [TJ since there do exist legit¬ 
imate black hole solutions here. For example, in the 
case where a is negative with k = — 1, for small p we 
see discontinuities in the Gibbs free energy which could 
be studied further to perhaps uncover new critical be¬ 
haviour. However, because this region exists beyond the 
pressure wedge of Figure[I] we cannot be certain that the 
branch we are considering has valid AdS asymptotics; it 
may be that this analysis jumps to branches for which 
f(r) terminates at finite r. Ensuring that the asymp¬ 
totics of the branches considered are well-behaved would 
require a non-trivial analysis of each branch, since the 
asymptotic structure will be a function of to, q, a and p. 
We shall not consider this case further. 


V. NEGATIVE p THERMODYNAMICS 


We now wish to consider the thermodynamic be¬ 
haviour of black holes which are solutions to the field 
equations ([9]) in the case of negative coupling. We wish to 
employ the same calculational machinery that was used 
in the previous section, but in order for this to be consis¬ 
tent we must substitute everywhere in our thermodynam¬ 
ics /i — > —p. This substitution allows us to explore the 
regions of negative coupling uses the same dimensionless 
parameters as before. The fundamental relationship for 
solutions of the field equations upon substitution p —> —p 
becomes, 


1 — ft T 



M k 3 


L 2 m 

r D-1 


q 2 L 2 

r 2(D-2) 


(42) 


where, as before, 


k = 


L 2 

r 2 


(F(r) - k) 


(43) 


with F(r) = r 2 g(r)/L 2 is defined so that the metric 
reads: — F(r)dt 2 + • • •. Similarly, any thermodynamic 
quantities can be obtained by substituting p —i —p in the 
aforementioned results. We can then employ the same 
dimensionless thermodynamic parameters as before 


i/a q l(D-2)(D-3) _ 

r h = vp' , q = 2 y--- » 

t = {D — 2)/i 1 / 4 T, a = X/y/p. 


(D- 3)/4 

5 


(44) 
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In terms of the quantities from eq. (44), the equation 
of state now reads, 


P = 


t (D-2)(D-3)k 2 akt (D — 2)(D — 5)fc 2 a 


v 

3k 2 t 


47 TV 2 V 3 

(D ~ 2)(-P — 7)k q 2 

47TV 6 ij2(D—2) 


47TU 4 


(45) 


and the dimensionless Gibbs free energy is given by 


9 = 


167T (u 4 
a(D — 8)v 
£>-4 


2 akv 2 — 3) 


3{D-2)kv D ~ 7 3(D — 2){D — 8)av D ~ 5 


4(D + 3)kv D ~ 3 


(D- 6) 


D—l 


607 rpv 


D -1 


(D-l)(D-6) 


— kv 


D+l 


(D — 6)(D — 4) D-6 

247rafcpu D+1 47rpv D+3 


2(D - 2)a 2 kv D ~ 3 
D - 4 


(D-l)(D-4) (D -1)(D- 2)J 


+ 


4(D — 3) (v 4 + 2 akv 2 — 3) ij 44 3 


(2D - 5)u 4 2(2D - 7)akv 2 3(2D - 9) 

D-2 + D-4 (D-6) 


(46) 


We once again restrict our attention to only those black 
holes which have positive entropy. In this case, the zeroes 
of the entropy are 


v± = 


W 


(D-2 


D-4 


—ka ± i a 2 + 


3(D — 4) 2 
(D — 2)(D — 6) 


(47) 

The positivity of entropy is straightforward for D > 7 
since the entropy will be positive provided v > u + for 
k = ±1. The situation is more complicated in D = 5 and 


is most simply captured by rewriting the bound from (47) 


as 


, 1 v 4 + 9 

ak > - - - 5 — 

0 v z 


(48) 


where v is the dimensionless specific volume defined ear¬ 
lier. 

We must also be aware of the asymptotic behaviour 
of the different branches of the solutions for this case. 


Substituting the dimensionless parameters from eq. (44) 
into eq. \42\ and taking the asymptotic limit, we find 
that the discriminant is given by 


A = 


167T 2 


4327 T z p 


2^4 


(D - 1) 2 (D - 2) 2 L (D — 2) 2 (D — l) 2 
167T o 72 tt 


-a 


(D — 1)(D — 2) (D — 1)(D — 2) 


a p 


+ (4 + a 2 )p 2 


(49) 


and is equal to zero for 

-9a-2a 3 ±2 (a 2 + 3) 3/2 . 

( 5 °) 

From this we see that p_ is strictly negative (and so 
does not yield asymptotically AdS behaviour) while on 



the other hand, p+ is strictly positive. A simple analysis 
reveals that A > 0 for p < p + and A < 0 for p > 
p + . From this we can conclude that in the region where 
p > p + , eq. (42) has one real branch and two complex 
conjugate branches, while in the region 0 < p < p+ all 
branches admit AdS asymptotics. Along the line p = 
p + all branches have AdS asymptotics and two branches 
coincide. 


As discussed in [3SJ US] , some AdS vacua are unstable 
in the sense that the graviton in these spacetimes may 
be a ghost. One can determine whether or not this is the 
case by examining the slope of the asymptotic form of 
eq. (|42|) 


h(h 


— 1 — kao + vttkA, + 


(51) 


where h'(jk oo) determines the sign of the kinetic term 
for the propagation of gravitons in the spacetime. The 
kinetic term has the correct sign when h' (/too) < 0 an d 
the ‘wrong’ sign when h! (kqo) > 0. First, we write h{k <*,) 
in terms of the dimensionless quantities, yielding 


k ( Tvqo ) 1 ^oo A 


and therefore, 


4a7 xpk\ 


167 j 2 p 2 k 3 


(D — 1)(D — 2) (D — 1) 2 (D — 2) 2 

(52) 


, 8 anpk os 48'K 2 p 2 K 2 00 

[hoo> (D-l)(D-2) (D-l) 2 (D-2) 2 ' 

(53) 

Expressing koo in terms of p and a by setting h(koo) = 0 
we find that h' (/Cqo) > 0 for p > p + for all branches, while 
there is one branch with h'(Koo) < 0 for p < p + . Hence 
the separation between stable and unstable regions coin¬ 
cides exactly with p + . The existence of AdS asymptotics 
and the stability of these regions is captured graphically 
in Figure [2] 
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FIG. 2. A plot showing the p — a parameter space in D = 5 
describing stability and existence of AdS vacua. The dark red 
curve is a plot of p +. For p > only one branch admits AdS 
asymptotics, while for 0 < p < p+ all branches admit AdS 
asymptotics. Also shown is the stability of the AdS vacua, 
for 0 < p < p+ there is one stable branch. For all values 
of D the plot is qualitatively the same, with the value of D 
providing an overall scale for the plot. 


Luckily for us, this stable AdS branch permits black 
hole solutions. The easiest way to see this is to solve 


, A , ll o L 2 m 

1 — k —- k 2 + —r = 


q 2 L 2 


L 2 


L 4 


„D -1 


r 2(D-2) 


where 


L 2 

K = g{r) -y k 


(54) 


(55) 


for g(r ) and note that a black hole solution would corre¬ 
spond to g{rh) = 0. We use the dimensionless parameters 
eq. (44) along with 


m = n ^ 3 ^ 4 m 


(56) 


as a rescaled mass parameter. A representative plot of 
g(r) is shown in Figure[3] Before moving on to thermody¬ 
namic considerations we note that the analysis performed 
here indicates that we are free to consider the thermo¬ 
dynamics of these black holes in the case a < 0, which 
corresponds to A, n < 0. This feature is captured in Fig¬ 
ure [2j where we see that there exist stable AdS vacua in 
this case. Furthermore, since there are no stable vacua 
for p > p + we regard this region of parameter space as 
unphysical. 



FIG. 3. A plot of the metric function g(r) demonstrating 
black hole solutions for the stable branch with p < p+. The 
red curve is for m = 25 and q = 2, the blue curve is for m = 30 
and q = 2. 


A. Thermodynamic Singularity 


A general feature that emerges in the case of negative 
coupling is the existence of thermodynamic singularities 
for both k = ±1. The k = 1 result is novel since in 
the case of Lovelock gravity, thermodynamic singularities 
were only found in the case k = — 1 [25]. Here we out¬ 
line some general considerations for the thermodynamic 
singularities that occur for negative coupling. 

As was discussed earlier, a thermodynamic singularity 
occurs when 


dp 

dt 



(57) 


Working with eq. ( |45| (first in the case where q = 0), 
we find that eq. (57) has the positive solution, 


= \j—ak ± \Ja 2 +)). 


(58) 


We immediately see that for k = ±1 only will be real, 
indicating only one thermodynamic singularity. We can 
define the thermodynamic singular point (p s ,u s ,f s ) by, 


, . v dp 

p s =p s (v s ,t s ) — 


= 0 ^ 
’ dv 


= 0. (59) 
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which yields 


= \J-ak + \/a 2 + 3 
D-2 


i s = 


127rafcuf + 27ruf — 307rw s 
47rq 2 1 


k(D — 3)w 4 + 2a(D — 5)t 


+ 3fc(L> - 7) + v2(D _ 5) 

(D-l)(D-2)k 


Ps = 


47 r {ak — Vet 2 + 3) {akVa 2 + 3 — a 2 — 3) 
ak( 2a 2 + 5) \/a 2 + 3 - 2a 4 - 8a 2 - 6) 


+ 


{ak — \/q! 2 + 3) 


D -2 ■ 


(60) 


Let us now make a few comments regarding the thermo¬ 
dynamic singularity. Studying the expression for p s in 
eq. (60) in the case k = — 1, we find 


Ps = P+ + 


{a + Vet 2 + 3) 


D -2 


(61) 



a 


and so in the hyperbolic case the thermodynamic singu¬ 
larity occurs at p s > p+ and is therefore excluded by our 
pressure and stability constraints. 

More interesting is the k = 1 case, where we have, 


Ps = P- + 


{—a + V oi 2 + 3) 


D-2 ■ 


(62) 


From eq. (50), p- < 0, so when q = 0 the thermody¬ 


namic singularity will occur for A > 0, and is therefore 
unphysical. However it will always be possible to choose 
q such that, for some range of a, p s will occur for physi¬ 
cally relevant pressures. This is illustrated in Figure Q. 
We will return to a deeper discussion of the thermody¬ 
namic singularity in the following sections. 


FIG. 4. A plot of p s for k = 1 for q = 0.8, 0.9,1.0 (blue, 
green, and red curves, respectively); the purple curve is p+. 
The thermodynamic singularity is physical for each value of 
q in the region where 0 < p s < p+. The plot is for D = 5 but 
the basic result is identical for D > 7 as well. 


1. Spherical 

We specialize first to the case of spherical horizons (i.e. 
k = +1) with q = 0. We find in this case that the number 
of possible critical points consistent with the constraints 
for a stable, asymptotically AdS vacuum breaks down as: 


B. Thermodynamics in D = 5 

In five dimensions the equation of state reads (for k = 

± 1 ), 


t 3 k 2 kat 3 1 


P = 


3fc _L q2 

+ -« (63) 


v 2nv 2 v 3 v 5 2 t rv e v e 

and the conditions for a critical point reduce to 

3 kv?; + 9k — 67rq 2 
c ttv c (v* + 6akv 2 — 15) ’ 

0 = — kv^ + 6v®a + ( 107 rq 2 — 60fc)t> 4 
+ (367rfcq 2 — 54 )au 2 — 307rq 2 + 45fc. (64) 


V 225 + 120VV) 

a < - zero possible critical points 


^225 + 120-/I5 

a = -—- single inflection point 


a G 


a > 


( V 225 + 120>/l5 V 282 + 102^17\ two possible 


\ 9 
V 282 + 102-\/l7 


6 I critical points 

one possible critical point. (65) 


At the point a = a/ 225 + 120-^15/9 the two inflection 
points ‘coalesce’ into a single inflection point. This inflec¬ 
tion point is not a physical critical point: the Gibbs en¬ 
ergy has a cusp at this combination of (p, t, v) , as shown 
in Figure ([5]), and therefore no phase transition occurs. 
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FIG. 5. Top : p—v plot for a = \/225 + 120\/l5/9 exhibiting 
an inflection point. Bottom,', g — t projection for p = p c , 1.1 p c 
and 0.9p c (red, blue, green, respectively). The Gibbs free 
energy has a cusp, indicating that the inflection point is not 
a critical point. 


More interesting is the situation that occurs when a £ 
(\/225 + 120\/l5/9, \/282 + 102-/L7/6). The behaviour 
of the Gibbs free energy and a p—v plot is shown in 
Figure [6j while Figure [7] shows a representative p — t 
plot. In this case there are two candidates for critical 
points. We find that there is a single physical critical 
point (characterized by the standard critical exponents), 
and the Gibbs free energy also shows a cusp. We find that 
for pressures and temperatures lower than and close to 
the physical critical point the Gibbs free energy displays 
classical “swallowtail” behaviour and the p—v plot shows 
typical van der Waals behaviour. 

The swallowtail is only present on the minimal branch 
of the Gibbs free energy for a small range of pressure 
and temperature, after which it is present on the up¬ 
per branch (see Figure [6]). After swapping to the upper 
branch, the tip of the swallowtail intersects the mini¬ 
mal branch, resulting in both zeroth order and first order 
phase transitions. This combination of zeroth and first 
order phase transitions leads to a reentrant phase tran¬ 
sition m, so named because after multiple (in this case 
two) phase transitions, the system returns to the same 
macroscopic type of state upon a monotonic change in 
the temperature (52] , 

To see this, first consider the lower left plot in Figure[6] 
For a temperature t < to the Gibbs energy is minimized 
for large black holes. Increasing the temperature, the sys¬ 
tem undergoes a zeroth order phase transition to a small 
black hole at t = t 0 as there must be a discontinuous 
jump to minimize the Gibbs energy. Further increasing 
the temperature, a first order phase transition occurs at 
t = ti, and the system is again a large black hole. In Fig¬ 
ure [7] we see the same effect - for pressures in the range 
p £~ (0.0057828,0.0059372), a large-small-large reen¬ 
trant phase transition is possible. The phase behaviour 
terminates at the point where the zeroth and first order 
phase transitions coincide in what is known as a “vir¬ 
tual triple point”. In all cases the critical behaviour is 
consistent with the maximal pressure constraint. 

In the region a > \J 282 + 102-\/l7/6 the single can¬ 
didate for a critical point is similar to the one just de¬ 
scribed, and the plots presented in the previous discus¬ 
sion are qualitatively identical to what is observed for 
these values of a. 
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FIG. 6. Reentrant phase transition: k = 1 case. Upper left-, p — v plot for a = 4, q = 0, k = 1. The two black curves 
correspond to the two possible critical points; only the upper black curve corresponds to a physical critical point. The blue curve 
is for t = 1.05f c , the red curve is for t = 0.96t c , and the green curve is for t = 0.83t c - Upper right : g — t projection for p = 0.95p c 
displaying classical swallowtail behaviour on the minimal branch, corresponding to a first order phase transition. Lower left : 
g — t projection for p = 0.92p c showing the swallowtail no longer present on the minimal branch. At to minimization of the 
Gibbs free energy requires a discontinuous jump, indicating a zeroth order phase transition. At t\ the swallowtail intersects 
the minimal branch and here is a first order phase transition. These two phase transitions combine to yield a reentrant phase 
transition. Lower right: g — t projection for p = 0.8p c . The swallowtail is on the non-minimal branch and there is no intersection 
with the minimal branch. There is no criticality in this region. 

We now turn our attention to the charged case. To 
begin with, we study the (q, a) parameter space to de¬ 
termine the number of critical points for particular com¬ 
binations of q and a. The results of this investigation 
are presented in Figure [8] There are no physical criti¬ 
cal points for a < 0. For a > 0, depending on the q, a 
combination there may be as many as two critical points. 


Figure [9] elucidates the thermodynamics in the case of 
two critical points. We present a representative analy¬ 
sis for the case q = 0.795. Here we observe an inter¬ 
esting feature for the critical pressure, namely that for 
a ~ 1.958731054 and a « 2.3834328996 two criti¬ 
cal pressures coincide, as can be seen in Figure [9] For 
a ~ 1.958731054 the two critical points coincide, while 
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FIG. 7. Reentrant phase transition: p-t plot for a = 4, k = 1. The figure on the right zooms in on the region of interest. 
The coexistence fine (black) starts at the virtual triple point and terminates at the physical critical point (green point). A 
zeroth-order phase transition (red) begins at the virtual triple point and continues until it meets with the curve marking the 
boundary for the ‘no black hole region’ (blue). We see that in this case there is a large-small-large reentrant phase transition. 


for a ~ 2.3834328996 only the critical pressures coin¬ 
cide, and the critical temperatures and volumes remain 
distinct. In each case, the critical point is governed by 
the standard mean field theory critical exponents. Hence 
these are not isolated critical points like that found in 


For a £ ~(1.95, 2.23) neither of the two possible crit¬ 
ical points are physical as they occur on a non-minimal 
branch of the Gibbs free energy. In this region we see an 
infinite coexistence line separating small and large black 
holes, as shown in Figure 10 For a > 2.23 we see the 
appearance of two separate physical critical points. A co¬ 
existence line beginning at the origin terminates at one 
of these critical points, while a coexistence line emanates 
from the second critical point, extending beyond physical 
pressures. As a is increased beyond ~ 2.23 these physical 
critical points occur for more widely separated pressures 
and eventually one of the critical points becomes unphys¬ 
ical - exceeding the maximal pressure constraint. 

For a small range of a £ ~(2.23, 2.38) we find a triple 
point. This feature is presented for a = 2.33 in Figure [TT| 
Referring to the p — v plot, we see the ‘double van der 
Waals’ oscillation that is characteristic of a triple point. 
In the p — t plane we see two coexistence lines for first 
order phase transitions that come together at the triple 
point. Near the triple point we see the presence of a 


small/intermediate/large black hole phase transition that 
is reminiscent of the solid/liquid/gas phase transition. 
Considering the Gibbs free energy, we see in Figure [yj 
that for pressures slightly above the triple point there 
are two first order phase transitions. As the pressure is 
decreased, these first order phase transitions come closer 
together until they finally merge at the triple point, as is 
shown in Figure [TTji. 


As discussed earlier, the k = 1 thermodynamic sin¬ 
gularity can occur for physical pressures when charge is 
included. One may wonder if the thermodynamic singu¬ 
larity ever coincides with a critical point. Simple calcu¬ 
lations show that provided 

1 / 36 — 2-7r 2 q 4 — 127rq 2 + 27rq 3 y / 7r(7rq 2 + 12) 

= 2 V 27rq 2 - 3 

( 66 ) 

then the thermodynamic singularity occurs at a potential 
critical point; an analysis of its features indicates that it 
is an isolated critical point. However we find that the 
maximal pressure constraint is always violated for the 
value of a at which the thermodynamic singularity co¬ 
incides with the critical point. Hence this this critical 
point is always unphysical. 
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a 



a 


FIG. 9. A plot of p c vs. a for D — 5, k = 1 and q = 0.795. 
We see that at the intersection of the red and blue curves 
(a « 2.3834328996) the two critical pressures coincide. The 
purple curve is the maximum pressure constraint. 



FIG. 8. Critical points in (a, q) parameter space for 

D = 5 and k = 1. Top: Possible critical points with 
( p c ,Vc,t c ) satisfying pressure and entropy constraints. Gray 
dots indicate zero critical points, blue dots indicate a single 
critical point, and red dots indicate two critical points. Bot¬ 
tom: A close up view of the top plot focusing on 0.7 < q < 0.9 
showing more detail. 


When the thermodynamic singularity is physical it ap¬ 
pears as a point at which all of the isotherms intersect at 
a particular volume. This feature is highlighted in Fig¬ 
ure 12 for a = 2.4 and q = 0.795. Evaluating eqs. (60) 


for the these parameter values we find that, 


p s « 0.02835105681, t s « 0.1267089635, 

v s ^ 0.7481508651. (67) 

Expanding the Gibbs free energy and temperature about 
the thermodynamic singular point we find that, 

QT) 

g(v s + dv,p s + dp) = A- —b B + 0(dp, dv) 
dv 

t(v s + dv,p s + dp) = C -j- + t s + 0(dp , dv) (68) 
dv 

where A, B and C are numerical constants. We observe 
that the Gibbs free energy expressed as a function of v 
and p does indeed diverge at the thermodynamic singu¬ 
larity. However, as noted in [2S], the Gibbs free energy 
is more naturally expressed as a function of t and p. In 
terms of these natural variables, the expansion of the 
Gibbs free energy near the thermodynamic singularity 
is, 

g(t s + dt,p s + dp) = —dt + B + (D(dp, dv) (69) 

o 

which is finite and smooth. We therefore conclude that 
the thermodynamics are well-behaved at the thermody¬ 
namic singular point. We note that when the thermody¬ 
namic singularity does not occur at a potential critical 
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FIG. 10. Charged case: p — t plots for q = 0.795, k = 1, D = 5. Left: p — t plot for a = 2 showing the absence of physical 
critical points and a coexistence line indicating a first order phase transition between small and large black holes. Center, p — t 
plot for a = 2.43 here we see two physical critical points (green dots). The lower coexistence line displays standard van der 
Waals behavior: beginning at (0,0) and terminating at a critical point. The upper coexistence line displays ‘reverse van der 
Waals’ behaviour: beginning at a critical point and continuing on indefinitely. Right: p — t plot for a = 2.5 showing only a 
single physical critical point and van der Waals-type behaviour. 


point (which is the case for the physical critical points 
considered here), the system is governed by the standard 
mean field theory critical exponents given by eq. (41). 


2. Hyperbolic 

In the case k = — 1 we find that, regardless of the value 
of q, for both positive and negative a there is only one 
possible critical point (positive v c ,t c and p c ). However, 
in all cases this critical point is unphysical, since we find 
p c exceeds the maximum pressure constraints. This be¬ 
ing the case, we still observe interesting thermodynamic 
phenomena. 


For a < 0, positivity of entropy is automatically satis¬ 
fied. Here we observe the presence of a large/small/large 
reentrant phase transition, as can be seen in Figure [13] 
which shows a representative p—t plot for a = —5 and 
q = 0. We find that the actual critical point is unphysical: 
it occurs at a pressure exceeding the maximum pressure 
constraint. However, the coexistence line that begins at 
this critical point dips below the maximum pressure con¬ 
straint, allowing for physical first order phase transitions. 
Since the coexistence line forms a “U” shape, it is possible 
for two successive first order phase transitions to occur, 
resulting in a reentrant phase transition. Note that a 
qualitatively similar p — t plot is found for all a < 0 and 
for q > 0. 

When considering a > 0, the positivity of entropy con¬ 
dition is no longer trivial and must be enforced by ensur¬ 
ing that eq. (481 is satisfied for k = —1. For a < ~ 1.000 


we observe features qualitatively similar to the a < 0 
case, and Figure [13] suffices to display the physics at play. 
However, as a becomes larger (still regarding black holes 
with negative entropy as unphysical), discontinuities are 
introduced into the Gibbs free energy. These disconti¬ 
nuities give rise to zeroth order phase transitions which 
lead to multiple reentrant phase transitions. This feature 
is illustrated in Figure pT| for q = 0 and a = 3. 

Referring to Figure |T4])a) we see how the structure 
of the Gibbs free energy gives rise to a multiple reen¬ 
trant phase transition. Due to the positivity of entropy 
constraints, there are only two physical branches of the 
Gibbs free energy, and these terminate at finite temper¬ 
ature. Thus up to a certain temperature, there are no 
physical black holes. For smaller pressures, we observe 
a zeroth order phase transition as there is a discontinu¬ 
ous jump required to minimize the Gibbs free energy as 
the temperature increases. For larger pressures, the two 
branches of the Gibbs free energy intersect twice giving 
rise to two first order phase transitions. At a larger pres¬ 
sure (p ~ 0.0581) the zeroth order phase transition coin¬ 
cides with the first order phase transition that occurs at 
the lowest temperature. Above this pressure we observe 
only a single first order phase transition. 

Now considering Figure [l4|(c) we see that for a fixed 
pressure p in the region (po,p\) a monotonic variation of 
the temperature will result in a zeroth order small/large 
phase transition followed by a first order large/small 
phase transition and finally a first order small/large 
phase transition. Thus we have a two reentrant phase 
transitions (small/large/small and large/small/large). 
An identical result was found in [25] for ID 3 rd - 0 rder 
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FIG. 11. Triple point: q = 0.795, a = 2.33, k = 1, D = 5. Upper left: p — v plot at the triple point. Upper right: p — t plot 
for the triple point. Note the presence of two first order phase transition coexistence lines distinguishing small, intermediate 
and large black holes. At the triple point the coexistence lines meet. The lower coexistence line terminates at a critical point 
(green dot), while the upper coexistence line continues indefinitely, eventually becoming unphysical when it crosses the maximal 
pressure constraint (greater than the values of p shown on the plot). Lower left: g — t projection for p > p$ c showing the presence 
of two first order phase transitions, indicated by the black dotted line. Lower right: g — t projection for p = p 3 C . At the triple 
point the two first order phase transitions merge together. 

Lovelock gravity. This result is the first example of mul¬ 
tiple reentrant phase transitions for black holes in fewer 
than 7 dimensions. 

By further increasing a we find that the range of pres¬ 
sure for which multiple reentrant phase transitions occurs 
decreases until eventually at a ~ 4.916 we cease to ob¬ 
serve multiple reentrant phase transitions. For a > 4.916 
we observe a single large/small/large reentrant phase 


transition, as shown in Figure [15] Again, the zeroth order 
phase transition is a product of the fact that negative en¬ 
tropy black holes are regarded as unplrysical, and hence 
discontinuities are introduced into the Gibbs free energy. 

As a increases further, we observe thermodynamic be¬ 
haviour that is qualitatively similar to what we previ¬ 
ously described. However, as a becomes larger, the ze¬ 
roth order and first order phase transitions begin at larger 
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FIG. 12. A p — v plot for q = 0.795, k = 1, D = 5, a = 
2.4 showing multiple isotherms intersecting at a single point, 
indicating the presence of a thermodynamic singularity. 



FIG. 13. A representative p — t plot using a = — 5, q = 0, 
k = —1 showing the presence of a large/small/large reentrant 
phase transition. The horizontal red line corresponds to the 
maximal pressure constraint. The coexistence line originates 
at an unphysical critical point (green dot), but the coexistence 
line itself extends into the physical regime for a large range of 
temperatures and pressures. The temperature axis has been 
represented logarithmically to highlight the presence of the 
reentrant phase transition. 


values of the pressure. The result is that at a ~ 12 all 
phase transitions occur for pressures larger than the max¬ 
imal pressure constraint. 

In the case of nonzero q we observe no qualitatively 
different behaviour. The charge simply acts to shift the 
phase transitions in both a and p. 


C. Thermodynamics in D > 7 

Employing the positivity of entropy condition, we find 
that for k = ±1 there are no critical points in D = 
7, 8, 9,10 that have p c < p+ regardless of the value of 
q. Furthermore, we find there to be no interesting phase 
phenomena present. Accounting for the positivity of en¬ 
tropy, the Gibbs free energy displays two branches that 
meet at a cusp for all values of a and q. 

Based on the form of the equation of state we expect 
this situation to hold for all D > 7. Furthermore, since 
the thermodynamics in quasitopological gravity and 3 rd - 
order Lovelock gravity are identical for D > 7, we expect 
the same results will hold if this analysis was performed 
for 3 rd -order Lovelock gravity with negative coupling. 


VI. CONCLUSION 

Our investigation of the thermodynamic behaviour of 
black holes (coupled to a 17(1) gauge field) in cubic qua¬ 
sitopological gravity has yielded a number of interesting 
results. We found that their thermodynamic behaviour 
was quite similar in a number of respects to their counter¬ 
parts in cubic Lovelock gravity, but also found a number 
of striking differences. Throughout our study we have 
ensured that the parameters respect the physical crite¬ 
ria of positive black hole entropy and stable ghost-free 
vacua via the implementation of various thermodynamic 
constraints. 

As far as the similarities are concerned, we found that 
the equation of state in D > 7 is identical to the cu¬ 
bic Lovelock case, and so p > 0 black holes will exhibit 
phase behaviour similar to that already studied for this 
case [55]. This means black holes in cubic quasitopologi¬ 
cal gravity exhibit the same novel behaviour, particularly 
isolated critical points with non-standard critical expo¬ 
nents gp . 

However the equation of state for the quasitopologi¬ 
cal case extends to D = 5, unlike the cubic Lovelock 
case. Our investigation of this sector indicated that for 
spherical black holes there is a single critical point (with 
standard critical exponents) provided a > 1.9909. In the 
hyperbolic case no critical phenomena were found. 

The preceding considerations hold for p > 0. For 
p < 0 our study of the thermodynamics of 73 = 5 
black holes proved to be particularly fruitful. For 
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FIG. 14. Multiple reentrant phase transition: q = 0, a = 3, k = —1, D = 5. Left: A plot of the Gibbs free energy in 
the parameter range where multiple reentrant phase transitions are present. Center, p — t phase plot. The solid black line 
represents the coexistence line associated with a first order phase transition, the solid red line represents a zeroth-order phase 
transition. To the left of the dotted blue line there can exist no physical black holes (those with positive entropy). The maximal 
pressure constraint is shown as the dark red horizontal line. Right A zoomed in view of the center figure showing the presence 
of a small/large/small/large multiple reentrant phase transition for pressures between po and p\. 




FIG. 15. Reentrant phase transition a = 10, q = 0, k = —1: Left-. A g — t projection for p = 0.0219. Enforcing the positive 
entropy constraint leaves only two physical branches of the Gibbs free energy, both of which begin at non-zero temperature, 
indicated a ‘no black hole’ region. Right: A representative p — t plot for showing the presence of a large/small/large reentrant 
phase transition. The horizontal dark red line corresponds to the maximal pressure constraint. A first order phase transition 
(black line) and a zeroth order phase transition (red line) are present. The dotted blue line marks the region where the black 
holes have negative entropy (and hence are regarded as unphysical). 
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k = +1 spherical black holes we found examples of 
small/large/small black hole reentrant phase transitions 
in the uncharged case, and the presence of a triple point 
and a small/intermediate/large black hole phase transi¬ 
tion reminiscent of the solid/liquid/gas phase transition 
in the charged case. Although both of these features 
have been previously found in gravitational systems (first 
studied in the context of the Kerr-AdS solution in six di¬ 
mensions 1223 ), this work represents the first instance of 
these features observed in five dimensions. All of these 
were obtained for ghost-free vacua satisfying constraints 
of standard AdS asymptotics and positive entropy. While 
we found an example of an isolated critical point, it oc¬ 
curred within a ghosty AdS vacuum. 

For k = — 1 hyperbolic black holes we found examples 
of small/large/small black hole reentrant phase transi¬ 
tions as well as small/large/small/large black hole mul¬ 
tiple reentrant phase transitions. This latter case rep¬ 
resents the first time a multiple reentrant phase transi¬ 
tion has been observed in a five dimensional gravitational 
system, having previously only been found in seven and 
higher dimensions [58] . We found all instances of isolated 
critical points to occur for negative pressures. 

We found no new or interesting thermodynamic phe¬ 
nomena in D = 7, 8, 9 or 10 for negative quasitopological 
coupling. Based on the form of the equation of state, 
we expect this situation to hold for all D > 7 and for 
3rd-order Lovelock theories with negative coupling. 

Further thermodynamic territory remains to be 
mapped out in quasitopological gravity. Our analysis 
could be extended to negative a for /j, > 0. While this 
region of parameter space contains black hole solutions 
and appears to have some critical behaviour, the anal¬ 
ysis is complicated by the fact that one is not guaran¬ 
teed to have AdS asymptotics for all branches. Hence 
it could be that the thermodynamics corresponds to a 
branch of a compact spacetime. To keep track of the 
asymptotic structure would require a nontrivial analysis 
of the asymptotics for all branches as a function of the 
parameters q, m, and a. 


ACKNOWLEDGMENTS 

We thank David Kubiznak for helpful advice and 
Maple resources. This work was partially funded by 
NSERC through their Discovery, Vanier, and PGS pro¬ 
grammes. 


19 


[1] J. Creighton and R. B. Mann, “Quasilocal 
Thermodynamics of Dilaton Gravity coupled to Gauge 
Fields,” Phys.Rev. D 52 (1995) 4569-4587, 

gr-qc/9505007 

[2] M. M. Caldarelli, G. Cognola, and D. Klemm, 
“Thermodynamics of Kerr-Newman-AdS Black Holes 
and Conformal Field Theories,” Class. Quant. Grav. 17 
(2000) 399-420, hep-th/9908022 

[3] D. Kastor, S. Ray, and J. Traschen, “Enthalpy and the 
Mechanics of AdS Black Holes,” Class. Quant. Grav. 26 
(Apr., 2009) 195011, 0904.2765 

[4] B. P. Dolan, “The cosmological constant and the black 
hole equation of state,” Class. Quant. Grav. 28 (Aug., 
2011) 125020, 1008.5023 

[5] B. P. Dolan, “Compressibility of rotating black holes,” 
Phys. Rev. D 84 (Sept., 2011) 127503, 1109.0198 

[6] B. P. Dolan, “Pressure and volume in the first law of 
black hole thermodynamics,” Class. Quantum Grav. 28 
(June, 2011) 235017, 1106.6260 

[7] B. P. Dolan, “Where is the PdV term in the fist law of 
black hole thermodynamics?,” 1209.1272 

[8] M. Cvetic, G. W. Gibbons, D. Kubiznak, and C. N. 
Pope, “Black Hole Enthalpy and an Entropy Inequality 
for the Thermodynamic Volume,” Phys.Rev.D 84 (Dec., 
2010) 024037, 1012.2888 

[9] A. Larranaga and A. Cardenas, “Geometric 
Thermodynamics of Schwarzschild-AdS black hole with 
a Cosmological Constant as State Variable,” J. Korean 
Phys. Soc. 60 (Aug., 2011) 7.987-992, 1108.2205 

[10] A. Larranaga and S. Mojica, “Geometric 
Thermodynamics of Kerr-AdS black hole with a 
Cosmological Constant as State Variable,” The 
Abraham Zelmanov Journal 5 (Apr., 2012) 68-77, 
1204.3696 

[11] G. W. Gibbons, “What is the Shape of a Black Hole?,” 
1201.2340 

[12] D. Kubiznak and R. B. Mann, “P-V criticality of 
charged AdS black holes,” JHEP 033 (May, 2012) 1207, 
1205.0559 

[13] S. Gunasekaran, D. Kubiznak, and R. B. Mann, 
“Extended phase space thermodynamics for charged 
and rotating black holes and Born-Infeld vacuum 
polarization,” 1208.6251 

[14] A. Belhaj, M. Chabab, H. E. Moumni, and M. B. Sedra, 
“On Thermodynamics of AdS Black Holes in Arbitrary 
Dimensions,” 1210.4617 

[15] H. Lu, Y. Pang, C. N. Pope, and J. Vazquez-Poritz, 
“AdS and Lifshitz Black Holes in Conformal and 
Einstein-Weyl Gravities,” Phys.Rev.D 86 (Apr., 2012) 
044011, 1204.1062 

[16] A. Smailagic and E. Spallucci, “Thermodynamical 
phases of a regular SAdS black hole,” 1212.5044 

[17] S. H. Hendi and M. H. Vahidinia, “Extended phase 
space thermodynamics and P-V criticality of black holes 
with nonlinear source,” Phys. Rev. D 88 (Dec., 2012) 
084045, 1212.6128 

[18] R.-G. Cai, L.-M. Cao, L. Li, and R.-Q. Yang, “P-V 
criticality in the extended phase space of Gauss-Bonnet 
black holes in AdS space,” JHEP 1309 no. 005, (June, 
2013) , 1306.6233 


[19] D. Kubiznak and R. B. Mann, “Black Hole Chemistry,” 
1404.2126 

[20] A. Chamblin, R. Emparan, C. V. Johnson, and R. C. 
Myers, “Charged AdS Black Holes and Catastrophic 
Holography,” Phys.Rev. D 60 (1999) 064018, 

hep-th/9902170 

[21] A. Chamblin, R. Emparan, C. V. Johnson, and R. C. 
Myers, “Large N Phases, Gravitational Instantons and 
the Nuts and Bolts of AdS Holography,” Phys.Rev. D 
59 (1999) 064010, hep-th/9808177 

[22] M. Cvetic and S. S. Gubser, “Thermodynamic Stability 
and Phases of General Spinning Branes,” JHEP 9907 
no. 010, (1999) , hep-th/9903132 

[23] M. Cvetic and S. S. Gubser, “Phases of R-charged 
Black Holes, Spinning Branes and Strongly Coupled 
Gauge Theories,” JHEP 9904 no. 024, (1999) , 
hep-th/9902195 

[24] B. P. Dolan, “On the thermodynamic stability of 
rotating black holes in higher dimensions - a 
comparison of thermodynamic ensembles,” 1312.6810 

[25] B. P. Dolan, “The compressibility of rotating black 
holes in D-dimensions,” 1308.5403 

[26] B. P. Dolan, D. Kastor, D. Kubiznak, R. B. Mann, and 
J. Traschen, “Thermodynamic Volumes and 
Isoperimetric Inequalities for de Sitter Black Holes,” 
1301.5926 

[27] N. Altamirano, D. Kubiznak, and R. B. Mann, 
“Reentrant Phase Transitions in Rotating AdS Black 
Holes,” Phys. Rev. D 88 (June, 2013) 101502, 

1306.5756 

[28] A. M. Frassino, D. Kubiznak, R. B. Mann, and 

F. Simovic, “Multiple Reentrant Phase Transitions and 
Triple Points in Lovelock Thermodynamics,” 

1406.7015 

[29] D.-C. Zou, S.-J. Zhang, and B. Wang, “Critical 
behavior of Born-Infeld AdS black holes in the extended 
phase space thermodynamics,” Phys. Rev. D 89 (Nov., 
2013) 044002,1311.7299 

[30] D.-C. Zou, Y. Liu, and B. Wang, “Critical behavior of 
charged Gauss-Bonnet AdS black holes in the grand 
canonical ensemble,” 1404.5194 

[31] M.-S. Ma and R. Zhao, “Phase transition and entropy 
spectrum of BTZ black hole in three-dimensional 
gravity with torsion,” Phys. Rev. D 89 (Oct., 2013) 
044005, 1310.1491 

[32] M.-S. Ma, F. Liu, and R. Zhao, “Continuous phase 
transition and critical behaviors of 3D black hole with 
torsion,” Class. Quantum Grav. 31 (Mar., 2014) 
31095001, 1403.0449 

[33] S.-W. Wei and Y.-X. Liu, “Triple points and phase 
diagrams in the extended phase space of charged 
Gauss-Bonnet black holes in AdS space,” 1402.2837 

[34] J.-X. Mo and W.-B. Liu, “Ehrenfest scheme for P — V 
criticality of higher dimensional charged black holes, 
rotating black holes and Gauss-Bonnet AdS black 
holes,” Phys. Rev. D 89, (Apr., 2014) 084057, 

1404.3872 

[35] J.-X. Mo, X.-X. Zeng, G.-Q. Li, X. Jiang, and W.-B. 
Liu, “A unified phase transition picture of the charged 
topological black hole in Horava-Lifshitz gravity,” JHEP 
56 no. 10, (Apr., 2013) , 1404.2497 





20 


[36] J.-X. Mo and W.-B. Liu, “P-V criticality of topological 
black holes in Lovelock Born-Infeld gravity,” Eur. Phys. 
J. C 74 no. 4, (Apr, 2014) 

[37] L.-C. Zhang, M.-S. Ma, H.-H. Zhao, and R. Zhao, 
“Thermodynamics of phase transition in higher 
dimensional Reissner-Nordstrom-de Sitter black hole,” 

1403.2151 

[38] H.-S. Liu, H. Lu, and C. N. Pope, “Thermodynamics of 
Einstein-Proca AdS Black Holes,” 1402.5153 

[39] Y. Liu, D.-C. Zou, and B. Wang, “Signature of the Van 
der Waals like small-large charged AdS black hole phase 
transition in quasinormal modes,” 1405.2644 

[40] C. V. Johnson, “Holographic Heat Engines,” 1404.5982 

[41] A. Rajagopal, D. Kubiznak, and R. B. Mann, “Van der 
Waals black hole,” Phys.Lett. B737 (2014) 277-279 
arXiv:1408.1105 [gr-qc] 

[42] T. Delsate and R. Mann, “Van Der Waals Black Holes 
in d dimensions,” JHEP 1502 (2015) 070, 

arXiv:1411.7850 [gr-qc] 

[43] N. Altamirano, D. Kubiznak, R. B. Mann, and 

Z. Sherkatghanad, “Thermodynamics of rotating black 
holes and black rings: phase transitions and 
thermodynamic volume,” Galaxies 2 (2014) 89-159 
arXiv:1401.2586 [hep-th] 

[44] B. P. Dolan, A. Kostouki, D. Kubiznak, and R. B. 
Mann, “Isolated critical point from Lovelock gravity,” 
Class. Quant. Gray. 31 no. 24, (2014) 242001. 

arXiv:1407.4783 [hep-th] 


[45] R. C. Myers and B. Robinson, “Black Holes in 
Quasi-topological Gravity,” Journal of High Energy 
Physics 2010 (Mar., 2010) 1-35, 1003.5357 

10.1007/ JHEP08(2010)067. 

[46] R. C. Myers, M. F. Paulos, and A. Sinha, “Holographic 
studies of quasi-topological gravity,” Journal of High 
Energy Physics 2010 (Apr., 2010) 1-44, 1004.2055 

10.1007/ JHEP08(2010)035. 

[47] J. Oliva and S. Ray, “A new cubic theory of gravity in 
five dimensions: Black hole, Birkhoff’s theorem and 
C-function,” Class.Quant.Grav. 27 (Mar., 2010) 225002, 

1003.4773 

[48] D. Lovelock, “The Einstein Tensor and Its 
Generalizations,” Journal of Mathematical Physics 12 
no. 3, (1971) 498-501 

[49] W. G. Brenna and R. B. Mann, “Quasi-topological 
Reissner-Nordstrom Black Holes,” Phys. Rev. D 86 
(June, 2012) 064035 1206.4738 

[50] R. B. Mann, “Pair Production of Topological anti de 
Sitter Black Holes,” Class. Quant. Grav. 14 (1997) 
L109-L114, gr-qc/9607071 

[51] D. Kastor, S. Ray, and J. Traschen, “Smarr Formula 
and an Extended First Law for Lovelock Gravity,” 

Class. Quant. Grav. 27 (2010) 235014, arXiv: 1005.5053 
[hep-th] 

[52] C. Hudson, “The mutual solubility of nicotine in 
water,” Z. Phys. Chem. 47 (1904) 113. 


